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SECTION  I 


INTRODUCTION  AND  SUMMARY 


This  third  quarterly  report  summarizes  the  progress  iTiadc 
during  the  period  between  1 April  1977  to  30  June  1977  in  the  VELA  Network 
Evaluation  and  Automatic  Processing  Research  program  being  carried  out  by 
Texas  Instruments  Incorporated  at  the  Seismic  Data  Analysis  Center  (SDAC) 
in  Alexandria,  Virginia.  The  five  program  tasks  are: 

• Evaluate  the  performance  of  the  Iranian  Long  Period  Array 
(ILPA)  and  the  Seismic  Research  Observatories  (SRO) 

• Develop  advanced  methods  for  detecting  seismic  events  and 
evaluate  the  event  detection  capability  of  postulated  networks 
of  seismic  detectors 

• Extract  long-period  seismic  event  waveforms  given  the  loca- 
tion and  origin  time  of  an  event 

• Develop  and  demonstrate  function  processes  for  interactive 
signal  processing;  standardize  and  automate  the  processing 
functions;  also  demonstrate  their  feasibility  for  the  graphics 
processing  of  events 

• Compute  instrument  and  path  corrected  spectra  of  a network 
of  long-period  seismic  sensors  and  estimate  seismic  source 
parameters. 

The  Research  Objective  Plan  drawn  up  and  approved  for  performing  these 
tasks  is  approximately  on  schedule.  During  this  quarter  we  presented  ten 
papers  at  the  Spring  American  Geophysical  Union  meeting  and  one  at  a special 
symposium  of  the  Institute  of  Electrical  and  Electronics  Engineers.  A 
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demonstration  of  our  interactive  seismic  processing  system  was  provided  for 
the  VELA  Seismological  Center  and  others.  Associated  data  items  were  de- 
livered. 

For  both  the  ILFA  and  SRO  data  evaluations,  work  was  concen- 
trated in  the  following  areas: 

• Expand  the  data  base  using  the  National  Earthquake  Informa- 
tion Service  (NEIS)  lists  available 

• Maintain  data  preparation  software  which  is  primarily  the  con- 
version of  edit  and  beamform  programs  to  the  Terminal  Sup- 
port (TS)  system 

• Prepare  software  for  extended  noise  analysis 

• Perform  the  routine  signal  processing  and  analysis  needed  to 
evaluate  each  station  separately  and  to  evaluate  all  stations  as 
a network  designed  to  locate  and  measure  earthquakes  and  pre- 
sumed explosions. 

To  develop  advanced  methods  for  detecting  seismic  events,  four 
research  tasks  are  being  pursued.  These  are  as  follows: 


Perfect  and  test  methods  for  automatically  timing  the  arrival 
of  seismic  events 

Determine  the  feasibility  of  using  an  adaptive  beamformer  to 
detect  weak  signals  gener.ated  by  seismic  events 

Develop  an  automatic  method  to  associate  the  detections  of  an 
event  received  at  a sub- set  of  seismic  stations  and  to  predict 
the  arrival  times  of  event  phases  at  all  stations  of  the  network 

Improve  the  efficiency  and  accuracy  of  estimating  a seismic 
network’s  capability  to  detect  and  locate  events  of  a given 
magnitude  and  from  a given  location. 
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To  develop  event  waveform  extraction  techniques,  two  tasks 
are  being  performed  as  follows: 

• Extract  long-period  tjodywaves  from  edited  single  sensor  three 
con:ponent  data  and  from  array  data 

• Cascade  several  tTiethods  of  long- period  event  waveform  esti- 
mations to  oj>timi7.e  the  extraction  of  weak  signals. 

To  develop  interactive  graphics  signal  analysis  capability,  the 
following  three  tasks  are  being  performed: 

• Standardize  seismic  processing  functions  for  performing  dis- 
crimination analysis 

• Extend  the  Interactive  Seismic  Processing  System  (ISPS)  by 
developing  a seismic  programming  language  to  generate  dis- 
plays, plots,  and  tables  from  data  and  information  files  stored 
on  the  disk,  and  to  perform  simple  analysis  functions  on  this 
data 

• Develop  software  for  transfering  data  to  a mass  store, 

rhe  final  program  task  is  to  derive  accurate  transmission  path 
corrections  and  source  parameters  from  long-period  data.  This  is  broken 
down  into  three  tasks  as  follows: 

• Use  newly  developed  and  existing  software  for  spectral  analysis 
corrected  for  instrument  response,  path  attenuation,  and  dis- 
persion 

• Form  a suitable  data  base  with  long-period  waveforms  mea- 
sured at  a large  number  of  stations  providing  adequate  azi- 
muthal coverage  of  the  event 
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Compare  these  to  theoretical  spectra  and  radiation  patterns 
derived  from  the  best  fit  source  model  as  a means  of  estimat- 
ing the  source  parameters  associated  with  an  event. 
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SECTION  II 

ILPA  AND  SRO  EVALUATION 


A.  CURRENT  STATUS 
1.  SRO 

Eharing  the  last  quarter,  our  major  emphasis  has  been  on  sig- 
nal processing  and  analysis.  To  date,  a total  of  890  long-period  events  have 
been  processed  and  analyzed,  with  an  additional  2E5  events  processed  but  not 
yet  analyzed.  Processing  and  analysis  of  850  short- period  events  have  been 
completed,  with  an  additional  250  events  complete  at  most  stations. 

Measurement  of  short-period  RMS  noise  in  the  0.  5 - 4.  0 Hz 
passband  has  been  carried  out  along  with  the  short-period  signal  analysis. 
During  June,  we  also  began  measuring  the  maximum  one- second  noise  am- 
plitude. 

During  the  last  quarter  we  have  also  created  and  checked  out 
long-period  noise  analysis  programs  which  will  compute  RMS  noise  levels 
(with  and  without  corrections  for  instrument  response),  RMS  noise  level 
trends,  maximum  25  second  noise  amplitudes,  and  noise  spectra. 

We  have  held  off  processing  long-period  noise  samples  until 
a new  unified  SRO  processing  program,  which  combines  our  existing  edit, 
data  preparation,  and  plot  programs,  was  completely  checked  out.  This  pro- 
gram will  reduce  the  amount  of  computer  time  we  will  use.  It  is  now  ready 
for  use.  We  also  have  available  a list  of  the  time  gates  we  intend  to  use  as 
noise  samples. 

Preliminary  results  on  short-period  noise  measurements  in- 
dicate a very  strong  correlation  between  noise  levels  and  distance  from  the 
coast,  whereas  much  less  correlation  was  observed  for  the  long-period  data. 
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ILPA 


Two  major  areas  have  been  emphasized  during  the  last  quarter. 
In  the  area  of  signal  processing,  a total  of  537  events  from  the  time  jjeriod 
5/1/76  - 1 1 /12/76  have  been  processed  and  analyzed  for  event  detections  and 
data  quality.  Surface  wave  magnitude  estimates  are  currently  Ijeing  measured 
for  those  events  which  have  been  detected.  Estimation  of  signal  coherency  lias 
been  initiated.  The  program  to  accomplish  this  was  converted  to  TS  during 
the  last  quarter  and  six  test  events  processed  by  it. 

In  the  area  of  noise  analysis,  the  programs  to  study  propagating 
noise  have  been  converted  to  TS  and  tested  during  the  last  quarter.  Further- 
more, a total  of  146  noise  samples  covering  the  period  5/1  /76  - 3/12/77  have 
been  processed  and  visually  examined  for  the  presence  of  unreported  events. 
Approximately  40  of  these  samples  show  signals  and  hence  will  not  be  includ- 
ed in  the  noise  analysis. 

Preliminary  results  indicated  that  about  0.  1 increase  of  magni- 
tude capability  was  obtained  by  the  array  over  a single  site.  This  was  based 
on  data  collected  right  after  the  array  was  installed  and  was  greatly  influenced 
by  malfunctions.  Better  results  are  expected  when  the  array  operations  im- 
prove. 

B.  FUTURE  PLANS 

During  the  next  quarter,  work  will  be  concentrated  on  the  fol- 
lowing areas: 

• Routine  signal  processing  and  analysis 

• Noise  analysis 

• Report  preparation. 
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SECTION  III 


DETECTION  METHODS 


A.  CURRENT  STATUS 

The  short-period  (SP)  detection  and  timing  algorithm  is  near- 
ing completion.  Several  detection  and  timing  criteria  have  been  tried  on 
NORSAR  single  site  data.  The  abundance  of  secondary  SP  signal  arrivals 
makes  it  unfeasible  to  use  moving- window  phase  regression,  which  was  suc- 
cessfully used  on  long- period  (LP)  data.  SP  detection  and  timing  will,  in- 
stead, be  based  on  envelope  detection.  However,  the  instantaneous  phase 
and  frequency  functions  may  serve  as  a detection  confirination. 

For  the  adaptive  beamforming  task,  the  evaluation  of  the  new 
ABP'  was  started  this  quarter  using  the  recorded  short-period  array  data  from 
the  Korean  Seismic  Research  Station  (KSRS),  First,  a strong  event  from 
Kamchatka  was  used.  Second,  a signal  from  a presumed  explosion  in  eastern 
Kazakh  was  studied  in  detail.  The  ABF  performance  with  this  event  was  in- 
vestigated in  different  passbands  and  with  various  convergence  rates.  Third, 
an  ensemble  of  thirteen  weak  events,  which  was  studied  under  Project  6705, 
was  also  reprocessed.  The  new  ABF  performance  appeared  to  be  consistent 
in  both  high  and  low  amplitude  signals.  The  preliminary  conclusion  is  that  the 
detection  capability  of  the  new  ABF  processor  is  equivalent  to  increasing  the 
number  of  elements  in  a conventional  beamforming  array  by  a factor  of  three 
or  four.  It  should  increase  the  detection  capability  of  arrays  by  about  0.  3 
units.  We  need  to  examine  more  data  to  confirm  this  preliminary  result. 

Under  the  Detection  Association  Processor  (DAP)  task,  depth 
dependence  was  put  into  the  automatic  event  location  algorithm.  Functions 
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were  constructed  for  more  accurately  generating  the  travel  time  and  ray 
parameter.  The  earth  model  to  be  used  to  generate  simulated  event  detec- 
tion bulletins  was  modified  with  the  above  improved  functions  to  more  ac- 
curately reflect  the  effect  of  depth  on  the  simulated  observations  of  seismic 
waves. 

R.  FUTURE  PLANS 

The  following  developments  are  expected  in  the  next  quarter: 

• Complete  evaluation  of  a detector  for  timing  of  the  onset  of 
short-period  P-waves 

• Complete  evaluation  of  the  new  ABF  algorithms 

• Develop  improved  detection  association  processor  (DAP)  soft- 
ware and  test  modifications  with  simulated  network  detection 
bulletins 

• Complete  development  and  demonstrate  new  network  capability 
evaluation  software 

• Prepare  reports. 
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SECTION  IV 


EXTRACTION  OE  EVENT  WAVEFORMS 
A.  CURRENT  STATUS 

An  adaptive  filter  was  developed  to  extract  long-period  body- 
wave  phases  based  on  their  expected  polarization.  Tests  of  the  processor  in- 
dicated the  need  to  modify  the  polarization  criteria  used  for  detecting  shear 
waves.  Tests  were  performed  using  the  Wiener  filter  cascaded  with  the  three 
component  adaptive  filter  (TCA)  to  extract  weak  surface  waves.  Preliminary 
results  indicated  approximately  0.  4 magnitude  improvement  in  detection  capa- 
bility of  three  component  surface  W'aves  can  be  expected  by  cascading  the 
Wiener  filter  and  the  TCA  filter.  Processing  of  events  is  continuing  in  order 
to  evaluate  improvements  in  the  detection  threshold. 

H.  FUTURE  PTANS 

In  the  next  quarter,  the  modified  processor  for  the  extraction 
of  bodywaves  will  be  evaluated.  The  processing  of  surface  waves,  using  the 
cascaded  Wiener-  TCA  filter,  will  be  completed  during  the  next  quarter. 
Preparation  of  the  report  will  be  continued. 
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SECTION  V 


INTERACTIVE  SEISMIC  PROCESSING 
A.  CURRENT  STATUS 

IXiring  the  past  quarter,  the  Extended  Interactive  SeisiTiic  Pro- 
cessing System  (ISPSE)  was  demonstrated  and  delivered  to  AFTAC  along  with 
system  and  user  documentation.  On  delivery,  ISPSE  featured  the  program- 
mable mode  of  operation  developed  during  the  first  half  of  the  contract  period 
together  with  a standard  processing  function  for  baseline  discrimination.  Also 
work  on  design  and  implementation  of  the  Interactive  Seismic  Programming 
Language  (ISPL)  (the  additional  feature  required  to  coniplete  ISPSE)  was  be- 
gun, ISPl,  consists  of  the  Interactive  Seismic  Programming  Langviage  Com- 
piler (ISPLC)  and  the  Interactive  Seismic  Programming  I.anguage  Interpreter 
(ISPLI),  ISPLC  in  turn  has  two  main  parts,  PASSI  for  editing  and  syntax  a- 
nalysis  whicli  has  recently  been  completed  (i,  e,  , all  associated  routines  de- 
signed, coded,  and  debugged),  and  PASS2  for  semantic  analysis  and  tnatrix 
generation  on  which  work  is  proceeding.  Also  a modification  to  the  Cepstrum 
method  is  b 'ing  incorporated  into  the  Short- Period  Earthquake/Explosion  Dis- 
criminator (SPEED)  processing  module.  This  module,  together  witli  user 
functions  compiled  via  ISPL,  will  permit  the  definition  of  a more  comprehen- 
sive standard  function  for  discrimination  processing  by  the  end  of  the  contract 
period. 

Under  the  KSRS  data  transfer  task,  the  long-period  and  short- 
period  reformat  programs  have  been  coded.  We  have  received  the  necessary 
information  for  constructing  the  KSRS  transfer  from  Ceotecii  and  VELA  Seis- 
mological  Center,  As  a preliminary  test,  some  LP  data  has  been  success- 
fully transferred  to  the  mass  data  store. 
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' B.  FUTURE  PLANS 

Currently  work  is  focused  on  the  design,  coding,  and  checkout 
of  routines  associated  with  the  PASS2  portion  of  ISPLC.  When  this  task  is 
complete,  similar  efforts  will  be  directed  toward  the  ISPLI  module.  Also, 
we  are  imbeddimj  ISPI,  in  the  ISPSE  system,  maintaining  SPEED,  and  de- 
signing high  level  statements  for  ISPL. 

The  data  structures  required  to  integrate  ISPL  into  ISPSE  are 
being  designed  and  constructed.  The  'procedure'  for  editing  and  interpreting 
routines  is  being  modified  to  accept  the  user  defined  functions. 

ISPSE  is  being  modified  to  give  the  seismologist  the  capability 
to  measure  the  maximum  amplitude  of  a signal  trace  for  a specified  frequency. 
Defaults  1.0  and  25.0  second  periods  for  short-  and  long-period  waveforms, 
respectively,  will  be  set  up.  The  analyst  may  specify  any  period  for  the  mag- 
nitude measurement. 

The  following  developments  are  expected  to  be  completed  dur- 
ing the  fourth  quarter  under  the  KSRS  transfer  task: 

• Transmit  a sample  of  raw  data  reformated  into  II  events  to 
the  mass  store  SP  file 

• C instruct  the  written  operational  procedure  for  the  operator 
for  long-period,  short- period,  and  coarse  status  transmission 

• Retrieve  and  analyze  the  transmitted  test  data  and  sulimit  the 
result  for  approval 

• Document  the  K.SRS  reformat  and  transmission  programs  and 
document  the  ILPA  programs. 
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Problems 


After  the  last  update  to  ISFSE,  the  ld)P  system  software  gen- 
erated a task  that  could  not  be  installed.  Currently  all  effort  is  being  applied 
to  determine  the  nature  of  this  problem,  A solution  to  this  problem  is  of 
paramount  importance  before  other  efforts  on  ISPSE  can  be  continued. 
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SECTION  VI 

SOURCE  PARAMETERS  FROM  LONG- PERIOD 
SURFACE  WAVE  DATA 

A.  CURRENT  STATUS 

1,  Softwiire  Development  and  Data  Prei).i  ration 

IXiring  the  early  part  of  this  qiiarter,  the  last  piece  of  com- 
puter program  needed  in  this  task  was  written  in  order  to  correct  the  observed 
surface  wave  spectra  from  the  Special  Data  Collection  Syste*m  (SDCS)  stations. 
Also  the  processing  of  the  edited  event  signals  using  the  PDP-15  comjiuter  was 
finally  comyjleted  for  all  selected  NTS  (Nevada  Test  Site),  EKZ  (eastern  Kaz- 
akh), and  PNE2  (peaceful  nuclear  explosion)  events. 

Z.  Analysis 

During  the  first  and  seconfi  quarter,  the  analysis  part  of  this 
task  was  done  whenever  the  yirocessed  data  were  availalde.  After  the  com- 
pletion of  the  data  preparation  during  the  eiirly  part  of  this  quarter,  analysis 

work  has  been  progressing  smoothly.  The  analysis  work  consists  of  three 
I 

main  areas:  1)  examination  of  the  observed  surface  waves,  (Z)  estimation 

of  the  surface  wave  travel  path  attenuation,  and  (3)  estimation  of  the  seismic 
source  parameters. 

The  surface  wave  group  velocities  for  the  travel  paths  from 
the  NTS,  EK/,,  and  PNE  to  all  available  observation  stations  were  examined 
I using  data  from  all  of  the  selected  NTS,  EK/.,  and  PNE  events,  F3y  com- 

’ bining  many  events,  surface  wave  group  velocities  for  each  travel  path  were 

\ obtained. 
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The  ol)Scrved  surface  wave  amplitude  spectra  of  the  selected 
NTS,  KKZ,  and  PNE  events  were  examined  for  the  peiiod  ranye  between  10 
and  50  seconds.  The  amplitude  spectral  ratios  of  the  Love  wave  to  Kayleigh 
wave  (LQ/I.R)  were  calculated  for  each  event  at  all  available  observation  sta- 
tions. These  ratios  were  compared  amon;,  the  NTS,  KK/,  and  PNF.  events 
for  any  indications  of  ^eoloyical  and  sou rce - orie nted  dependence.  .Mso,  the 
relative  excitations  of  the  Love  waves  for  a suite  of  NTS  events  were  examin- 
ed in  terms  of  the  ratios  of  the  observed  Love  wave  amplitude  spectra  of  a 
particular  NTS  event  to  a reference  NTS  event. 

Preliminary  results  indicate  that  LQ/LR  ranges  from  zero  to 
ten  for  NTS  events. 

The  average  surface  wave  attenuation  coefficients  along  tiie 
travel  paths  between  the  NTS  and  the  observation  stations  and  those  along  the 
travel  paths  between  the  EKZ  and  the  region- station  attenuations  were  esti- 
mated by  the  isotropic  - source  method  using  all  available  events.  The  travel 
paths  between  the  NTS  and  the  observation  stations  were  logically  divided  into 
two  groups:  oceanic  path  and  continental  path.  In  addition,  wherever  applic- 
able, the  two- station  method  was  employed  to  calculate  the  surface  wave  at- 
tenuation coefficients  for  the  travel  path  between  two  observation  stations. 
These  result"  were  compared  with  the  existing  estimates  of  surface  wave  at- 
tenuation. 

Preliminary  results  on  ocean  path  attenuation  at  NTS  are 

""4 

1. 8 X 10"  km  at  10  seconds  period  with  gradual  decrease  to  1 . 1 y 1 0 at 

50  seconds.  For  continental  paths  to  NTS  and  EKZ  the  attenuations  were 

from  4x10  at  1 0 seconds,  to  1 . 6 x 1 0 at  20  seconds,  with  gradual  de- 

.4 

crease  to  0.  8 x 10  at  50  seconds. 

The  seismic  source  parameter  estimations  for  all  selected 
NTS,  EKZ,  and  PNE  events  were  completed.  Theoretical  radiation  patterns 
based  on  the  estimated  source  parameters  were  generated  for  several  events 
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with  a large  number  of  observation  stations.  The  observed  values  were 
superimposed  to  make  comparisons  between  the  observations  and  the  theo- 
retical radiation  patterns. 

Preliminary  results  indicate  that  NTS  induced  earthquakes  are 
nearly  vertical  strike-slip  with  relative  strength  ranging  from  0.5  to  1.0.  An 
EKZ  induced  earthquake  was  dip- slip  with  a relative  strength  of  0,56, 

B.  FUTURE  PLANS 

The  analysis  work  needed  in  this  task  is  practically  done.  All 
these  results  will  be  organized  and  the  writing  of  the  first  draft  of  the  tech- 
nical report  for  this  task  will  be  started. 
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